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^ ■ Abstract 

^ ■ During several years the data on different channels pp {J^^ = ^) 

^ ' — > 3 mesons presented by Crystal Barrel Collaboration were suc- 

ly-^ . cessfully analyzed by extracting the leading amplitude singularities 

O ! — pole singularities — with the aim to obtain information about 

two-meson resonances. But these analyses do not take into account 

0> ' three-body final state interactions (FSI) in an explicitly correct way. 

^ ' This paper is devoted to the consideration of this problem. Here I 

demonstrate how the coupled three-body equations may be written 

Qh! for the 7r°7r°7r°, 'qn^n^, rirjn^, XiC7r° channels in the annihilation 

I at rest using three-body dispersion relation N/D-method. 



1 Introduction 

During several years Crystal Barrel Collaboration presents the high statistics data on 
three meson production from the pp annihilation at rest. These data were successfully 
analyzed (see, for example, [|l], H, |^, ^) with the aim to search for new meson resonances in 
the region 1000-1600 MeV using the K- matrix formalism or a simplified dispersion N/D- 
method. There is a strong expectation based on QCD as well as on lattice calculations 
that the lowest scalar glueball is located in this region. Thus, the identification of 
scalar resonances in the mass region 800-2000 MeV and their classification in gg-nonets 
should be done to trap the lightest scalar glueball. The quark/qluon structure of these 
resonances may be determined from the analysis of coupling constants of these states to 
pseudoscalar mesons 0. The question is: is the K-matrix approximation - or simplified 
N/D-method - sufficient for this purpose? 

In this paper the three-body dispersion N/D-method is presented, which is based on 
the two-body unitarity condition and allows one to take into account FSI of three mesons. 
The basic principles of this technique were developed in 1960s and they were applied to 
the calculation of FSI in the K-meson decay into three pions (see [|, and references 
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therein). Recently this technique was used for the calculation of the — > Stt decay ITT 



and 0-meson production in the pp annihilation at rest |T2|. 

This paper is organized as follows. In Section 2 the two-particle discontinuity of the 
decay amplitude is written out, and the integral equation is derived. It is shown how 
to take into account not only the S-wave binary interactions but also interactions with 
higher angular momentum (P- and D-wave). In Section 3 this technique is generalized 
for the resonance and non-resonance two-particle interactions. The complete three-body 
dispersion equation is written, that may be used for the analysis of various three-particle 
reactions. In Section 4 it is shown how the coupled three-body equations could be written 
for TT^Tr^Tr", r]n^n^, t^t^tt^, KKtt^ channels in the pp annihilation from the ^ state. 



2 Two-particle discontinuity of the decay amplitude 



Let us start with the consideration of the decay of a scalar particle with the mass M and 
momentum P into three scalar particles with masses mi, m2 and and momenta ki, 
k2 and ^3. I would like to show that the energy dependence of this transition amplitude 
may be derived from the unitarity and analyticity. The imaginary part of the amplitude 
An for the transition 1 — n is given by 

ImMn = Wi2nr6iP-f2h)T:,nMn', (1) 

^ n' i=l 

where r„/„ is the scattering matrix. Equation (|l]) may be simplified if the approximation is 
used where i) binary interactions of particles are taken into account, three-particle forces 
being neglected; ii) there are no transitions through intermediate states with more than 
three particles. It is also supposed that there are no transitions through intermediate 
states involving other particles; still, this is not a crucial restriction. Below it is shown 
how the production of new particles in the intermediate state may be taken into account. 
In this approximation one may write eq.(P as follows: 

ImM, = ^{2n)^6{P - k, - - k^)T;^M^. (2) 

Three-particle forces being neglected, the amplitude M3 for the three-particle produc- 
tion is a sum of four terms: i) the direct production amplitude A(si2, S13, S23), which is 
assumed, for the sake of simplicity, to be free of singularities; ii) the amplitude ai2, where 
the last interaction is that of particles 1 and 2; likewise, there are similar terms Hi) 013 and 
iiii) 023. I take into account FSI with different orbital momenta L, hence the amplitude 
ttij is written as 

a,, = Y.FLA^Mj). (3) 

L 

Here Sij is the two-particle invariant mass squared, Sij = {ki + kj)"^. The function 
defines the angular distribution of decaying particles. In the case of the S-wave binary 
interaction Fq = 1, but for the higher- wave interaction the amplitude depends on the 
angle between interacting and spectator particles as well. 
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2.1 S-wave interaction 



Let us begin with the simplest case of the S-wave pair interactions. The decaying amph- 
tude is given by 

M{Si2, Sl3, S23) = A(Si2, Si3, S23) + A°2(Si2) + A^^is^s) + AI.,{s23) . (4) 

The connected part of the matrix T33 is not known, so I cannot directly calculate the 
contribution from three-particle intermediate state into the imaginary part of M and 
write down corresponding dispersion integral. So, I will explore the two-particle unitarity 
condition to derive the integral equation for the amplitude A^j. This method was described 
in details in ref . . The idea of such an approach is that one should consider the case of 
small external mass M < nii +m2 +^123, when only the scattering reactions are physically 
possible. A simple relation can be written in this case, and then I perform analytical 
continuation of the final equation over the mass M back to the decay region. Let us stress 
that the obtained expression differs from the contribution from the disconnected part of 
T into imaginary part M: this result satisfies the three-body unitarity condition, in which 
three particles interact by pairs only ||13|| . 



So, I write down an ordinary unitarity condition for the scattering in the channel of 
particles 1 and 2. This means that (M + ms) < Su is assumed. The discontinuity of the 
amplitude is equal to 

disCi2M{si2,Si3,S23) = 



= \j rf$12(^l, k2) (a + A?2(Si2) + A\,{sr3) + <(S23)) A°l2(^12), (5) 

where (i$i2(fci, /C2) is the phase volume of particles 1 and 2 

^$12(^1, k2) = {27^f5\P -h- k2) ^-^^S{ml - kl)5{ml - kl). (6) 

^2^2 is the S-wave two-particle scattering amplitude, which can be written in the disper- 
sion N/D method as a series 



Al^2{s) = G^is)G^is) + G^is)Bois)G^is) + Gl^is)BUs)G^is) + ... 

G(^{s)G^{s) 



(7) 



where Gq and Gq are left and right vertex functions and B{s) is the dispersion represen- 
tation of the loop diagram: 

Bo{s) = r ds' G^is')p,,is')GM ^ 



where pij is the two-particle phase space 



Piji^) = tf~a/[s - (mi + m2)][s - (mi - 7722)] . (9) 
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Vertex functions contain left-hand singularities related to the t-channel exchange dia- 
grams, while B-function has singularities due to the elastic scattering. It is not specified 
from the consideration of the scattering amplitude Ag^g of ^q. (^) whether both vertices 
Gq and Gq have these singularities or only one of them. In the case of three-body decay 
the situation is quite opposite. On the first sheet the decay amplitude has only singulari- 
ties at Sij = {mi + rrij)^, which are associated with the elastic scattering in the subchannel 
of particles i and j. This means that the vertex Gq is analytical function. The simplest 
choice is 

G^ = 1 (10) 

Now return to eq. (^. Since the functions A do not have the two-particle threshold 
singularity and thus do not have a discontinuity, the Ihs of eq. (^ is 

disci2M{si2, si3, S23) = disci2A°2{si2) . (11) 

As was stressed in ref. |T^, only one rescattering of particles 1 and 2 can be considered 
in the final state and a full set of binary rescatterings can be taken into account multi- 
plying by (1 — Bq{si2))~^. Thus, the two-particle discontinuity in this special case of one 
rescattering is defined as 

disCi2 A%{Si2) = \j C^$12(A;i,A;2)(^A(si2,Si3,S23) + A°3(si3) + A°3(s23))g'^(si2). (12) 

It is convenient to perform the phase space integration in eq. (|12D in the center-of-mass 
system of particles 1 and 2. In this frame 

Si3 = ml + ml- 2kiok3o + 2213 \ h \ \ k^ \ , 
fcio = ^i^^^^l^ , \h\=^/H^i, (13) 

= -^12 + ^3 £ 1^1= 1/ k^n - 



30 — o , ) I '^3 |— V ft-30 ^ '"-31 



where 2:13 = cos9i3, and 61^ is the angle between particles 1 and 3 in the cms of particles 
1 and 2, and s = M^. The expression for S23 is obtained from eq. (|1^) by the replacement 
1^2. From eq. (|T2|) one has: 

dlSCi2 ^?2(Sl2) = (As(Si2) + (^?3(S13))0 + (^23(523))o) G'^(si2)pi2 (^12) , (14) 

where the following notation is used: 

(^°3(^.3))o= f ^ ^^A%is.s). (15) 
Xs is the S-wave projection of A: 

As(si2);= y"^^^ A(si2,Si3,S23)- (16) 

Analytical continuation over external mass M from the scattering to the decay region 
allows one to define correctly the rules of integration over z. This integration should be 
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carried out along the contour Cj(si2), whose position at different S12 is described in detail 
in ref.[^. Here I would like to note that only at small S12, 

(mi + 7712) < S12 < \ (mi + m2 - mj - m^ma, (17) 

rrii + 1713 ""^i + ms 

it coincides with the phase space integration contour 

- 1 < ^^3 < 1, (18) 

and it contains an additional piece at larger S12. 

Equation ([I^ allows us to write down the dispersion integral for the amplitude with 
one pair rescattering in the final state: 

/(mi+m2)2 TT s'12 — S12 

Here I exclude A5 from the dispersion integral, but it is also possible to include it: in both 
cases the unitarity is satisfied but with different behaviour of the amplitude at the infinity, 
which cannot be defined by the unitarity and analyticity only. After the transition from 
one rescattering to the full set of binary interactions in the final state one has 



AUsu) = Xs{su)B{s,,) + r ^^^#^^^^(^3(4)0 + (A°3(4)o). (19) 

J(mi+m2)^ 71 S12 — S12 



.0 / N Ag(gi2) ^0(^12) , 

^^^^'"^ - 1 - 5o(.i.) + 
1 r ds',, /^-(^'-)^o (^'-) ((^0^(4))^ + «(4))o). (20) 



1 - i?o(Sl2) J(mi+m2)2 TT ^'12 " S12 

Let us now check that the extraction of the final state interaction does not violate the 
unitarity condition (|). To calculate the Ihs of eq. (|^), the equation (pOD is rewritten as 
follows: 

AUsu) = |^^~ff"\2 (^'(^-) + Asu)) , (21) 

1 1 - -80(^12) r V J 



where 



b\su) = Xs{si2) Boisu), (22) 
J{su) = r , ^ ^^^(fi^)^o(^i^) ((A°3(4))o + «(43))o). (23) 

J(mi+m2y 71 S12 — S12 



I J- — -Do (-512) I 



Thus, 



+ (1 - 5o(si2))/m b\su) + Im Bo{su)J{su) + (1 - Bo{si2))disc J(si2) j . (24) 
Taking into account that 

disc So(si2) = Im Bq^su) = ^12(512)^0(512), 

disc b^isu) = puisu) Gl^{si2)Xs, (25) 

disc J(Si2) = Pl2(Sl2)Go (Sl2)(^?3 + ^23)0, 

I get the rhs of eq. (j|), hence, the unitarity condition is fulfilled. 
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2.2 P-wave interaction in the final state 



First, determine the structure of the amphtude, where particles 1 and 2 interact in the P- 
wave and the particle 3 is spectator. The easiest way to do this is to transform the decay 
amplitude into scattering amplitude. To perform this transformation the antiparticle 3 
with the momentum (— ^3) is considered instead of particle 3 itself. The assumed form of 
the amphtude is Oi^Qif^A\2{si2)- Operator Qi describes the P-wave angular distribution 
of particles 1 and 2 in the final state, it is defined as relative momentum of these particles 

2 2 

QiM = fciM - - — {h + k^)^. (26) 

■S12 

Operator Oi should be constructed as relative momentum of initial state and antiparticle 
3. Taking into account that Qi^ + A;2)^ = 0, Oi is defined as 

Oi^ = h^,. (27) 

As is seen from eq. (^, 

Fi = Oi^Qi^, (28) 

and it is easy to find out that Fi is proportional to 2:13 in the cms of particles 1 and 2. 

Hereafter the procedure of the above section is used to write down the two-particle 
unitarity condition and the integral equation for A^. The P-wave two-particle scattering 
amplitude of particles 1 and 2 in the N/D method can be written as 

Al^,{su) = Qi, . ^'^'f , Qi„ (29) 

i — -Di(Si2j 

where the B-function is equal to 

ds' Gf{s')pijis'){Qi^Q,^) 



J(mi+mj)'^ TT S — Sij 

and (...) means the averaging over space angle. 

As the first step, one should consider triangle diagram with one P-wave rescattering of 
particles 1 and 2. The discontinuity of this diagram is equal to: 

discuAl^i^su) = ^ J d^i2{ki,k2)ai3{si3,zi3)Qi^G^{si2) = 

= G'f(Si2)pi2(Si2) y ^Ql/.ai3(si3,^13)- (31) 



The integration over space angle is performed in the cm. frame of particles 1 and 2. In 
this frame Qi turns into ki2. The z-axis being directed along ks, the components of ki2 
are equal to 

ki2x = ki2 sin6'i3 cos(/), ki2y = k^ sin 6^13 sin^, k^z = ku cos ^13 



ki2 = \ — [■S12 - {rrii + m2)2][si2 - (mi - ma)^], (32) 
V S12 
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where is azimuthal angle. The integration over (f) in eq. (|3T|) keeps the components of 
Qi^ with fi = z only: 

Qix • ■ • = / -j—Qiy ■ ■ ■ = 0, 



47r J 47r 



i2 • • • = ki2 J zi3... (33) 



Let us introduce the vector with only z component unequal to zero in the cms: 

ki = h,-p,^^, (34) 



where p = ki + k2- Then it follows from eq. ([ 

disc = kLpi2{si2)G[{si2)—ri^= I ^^213013(513,2:13), (35) 

where 

_ ^^±Y = [(^-^3)^-gi2][(M + m3)^-Si2] ^3g^ 

^ 4Si2 

The invariant part of the P-wave amplitude can be written as a dispersion integral over the 



energy squared of particles 1 and 2 in the intermediate state, while k^ in eq. ([35| ) defines 
the operator structure of the P-wave amplitude. Thus, the P-wave triangle diagram is 
equal to 

,± r Pl2{s'i2)Gi{s\2) I ^, \ (^^. 

^J(mi+m2)^ IT S12 - S12 

where 

/ \ ^12 f dz . . . . 

(013)1 = / / ^2; ai3(si3,z). (38) 

^-{kif Jcz{s,2) 2 

Equation (137|) should be multiplied by the operator Qi which describes angular distribu- 
tion of particles 1 and 2 in the final state. Using /c^ Qi^ = Oi^ Qi^, one may conclude 
that the operator part of 012 is correctly reconstructed. To take into account binary rescat- 
tering in the final state it is necessary to multiply eq. ( P7D by the factor (1 — -81(512))^^. 

The same steps should be done, if 013 is replaced by the direct production term A(si2, S13, S23) 
Still, the energy dependence of A is not usually known, so use a simpler assumption is 
used. Let us introduce the constants which define the direct production amplitude of 
particles i and j in the P-wave. Then the following integral equation for the amplitude 
can be written: 



. w X Af2 Si(si2j 
Al2{Sl2) 



1 (^s'i2pl2(s'i2)Gf'(s'i2 



1 - i?l(Si2) 7(mi+m2)2 TT s'12 - S12 



5l(Sl2) 

■^\(ai3)i + (023)1). (39) 
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2.3 D-wave interaction in the final state 

Likewise the case of D-wave interaction in the final state may be investigated. The 
amphtude, where particles 1 and 2 have the last interaction in the D-wave, has a form 
O2 (52^12(5 12)- The operator Q2 is a traceless tensor of the range 2 constructed of relative 
momenta of particles 1 and 2: 

Q2fj.u = ki2^ki2u — -^I25'^i/, (40) 

where 

2 2 

ffl — Tfl 

ki2t^ = hf, - k2^, {ki + A;2)^, (41) 

■S12 

and 

9t = 9,. - (42) 
Another operator, O2, can be defined as 

02f,u = k-Sf^ksi,. (43) 

Below it will be proved that O2 can be defined in this way. Let us note here that the 
operator part of the D-wave decaying amplitude is proportional to 3zf^ — 1. 
Let us define the D-wave scattering amplitude of particles 1 and 2 in N/D representation 

as 

^^2(^12) = Q2,u T^^^ Q2,u. (44) 

The discontinuity of the triangle diagram with the D-wave rescattering of particles 1 and 
2 is equal to 

discuAl^isu) = ^ J d^i2{ki,k2)ai3{si3,Zi3)Q2f,uG2{si2) = 

= G2iSi2)pl2isi2) J ^Q2f.,.ais{si3,Zi3). (45) 

The integration over space angle is performed in the cm. frame of particles 1 and 2. 
Equation ( |32D is used, and only the following components of tensor kufikuu are not equal 
to zero: 

— K^UxK^Ux ■ ■■ ~ J ^ f^UyKUy ■■■ — J U " ^IsJ • • • 

, 7 7 2 f '^'^13 2 



kl2zkl2z ■■■ = k-^^J . . . (46) 



To perform the space integration in the invariant form, let us define the tensor (yf^^ as 

p2 {k^ 
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In the cm. frame of particles 1 and 2 the tensor gj^^^ has only two non-zero diagonal 
elements: 

g^^ ^diag{Q,-l,-l,Q), (48) 



Thus, eq. (^) can be rewritten as 

— ki2,iki2u . . . = -gfi^ I — (1 - 

k'u / % • • • (49) 



k k - n^^ h'^ [ ^^^^ n ^2 1 

kj^skus , 2 f 



{ki 



\2 



Equation ( ^91 ) should be placed into eq. (|i5| ) for the discontinuity of triangle diagrams. 
Taking into account that this expression must be multiplied by the external operator Q2, 
one can see that only the third term in eq. ( PP] ) contributes into eq. (^5]). Moreover, one 



could replace k^ by k^. Finally, as follows from eq. (^5]), 

disc = /C3/.A;3j.Pi2(S12)G2(Si2) (013)2, (50) 

where 

(^13)2 = jrr^ / ^ — (1-3^2) ai3(si3,2;). (51) 

1% ) "'C'3(S12) ^ 

The invariant part of the D-wave amplitude can be written as a dispersion integral and a 
full set of binary rescatterings in the final state are defined by the factor (1 — -B2(si2))~^. 
As a result, I come to the following integral equation for the amplitude Al2'- 

.2 / N _ Af2 ^2(^12) 

"^^^^'"^ - 1-52(S12) + 



ds'u Pi2(s'i2)G'^(s'i2) 



1 - i?2(Sl2) "'(mi+m2)2 TT s'12 - S12 



((ai3)2 + (023)2), (52) 



where the constant Af^ stands for the direct production amplitude of particles i and j in 
the D-wave. 



3 Dispersion equations with resonance and non-resonance 
production of particles 

In this Section more realistic case is considered, when there are resonance and non- 
resonance interactions between two particles. This situation happens, for example, in 
the O"*"^ wave of the pion-pion amplitude, where there is a non-resonance background and 
a set of resonances at the energies above 1 GeV. 

I start with the dispersion representation of the two-particle amplitude for this particular 
case. As the simplest case, consider the S-wave one-channel amplitude. The first resonance 
term of the amplitude can be written as 



(53) 
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where is a pure (non-physical) mass of resonance a, and the function g°'{s) describes 
its decay into two particles. Generally, these functions can depend on s, so the following 
form is suggested: 

g^is)=g''<j>is). (54) 
The second term of the amplitude with one virtual loop is equal to 

where b{s) is loop diagram with the cutoff function. The behaviour of (p{s) at large s (or 
at small distances r) is not known, and the simplest way to avoid these uncertainties is 
to assume that the contribution from r < tq is equal to zero. Then the cutoff function 
A(s) is defined as follows: 

A(.) = J d'r e^%(r - ro) / ^ e"^''U\s'), (56) 

and b{s) is equal to 

b{s) = tp{s) + P f— 4^ Ms')- (57) 

J IT S — S 

Summing up the terms with different number of loops, one obtains the following expression 
for the amplitude 

A= ^" , . (58) 



For non-resonance interaction of particles defined by vertex function f{s), eq. (|58|) should 
be rewritten as follows: 

A= ^ ^, (59) 



where 



bj{s) 



ip{s) + P l-4^Kf{s'), (60) 

J 71 S' — S 



and Af{s) is given by eq. (p6D with the replacement of (s) by f{s). It should be noted 
here that if and / are constant, then cutoff functions are equal to zero, and dispersion 
representation coincides with the K-matrix approach. 

Equation (pD| ) can be easily generalized for the case of the two-particle interaction with 
orbital momentum L: 

Al = r . Q.., (61) 



where b{s) and 6/(s) are given by eqs. (|57D and (^), with the cutoff functions A(s) and 
A/(s), correspondingly, which are equal to 

A(.) = j d'r e^%(r " ^o) / ^ e-^^V(s')(QL,,gL,,), (62) 
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A/(^) = / d'r e^%(r - ^o) / ^ e-^''^y{s'){Q,,,Q,,,) (63) 

The same formalism can be applied to the case of multichannel amplitude. The decay 
of resonance a into particles m and n is given by the function g!;^l{s) = gmn^rnnis) and 
the non-resonance transition from the channel with particles m and n to the channel with 
particles m' and n' is given by fmn;m'n'{s). Here the expression for multichannel amplitude 
is not given in an explicit form, the recipe of its construction can be found, for example. 



in [0]. I would like only to note that the denominator in eq. (0), which describes the 
rescattering of particles in multichannel case, has the matrix form (/ — B)~^, where / is 
a unit matrix and the B-matrix element is equal to 

Bmn;m'n' ^mn\m' n' {s\ / „ -|- Pmn\m'n'{s) , (64) 

J n s — s 

(3mn;m'n'{s) = ipmn{s) + P [ Pjipi — 1 Af^rnn;m'n'{s')- (66) 



Let us generalize the results of the previous sections and write down the dispersion 
equations for the three particle decay amplitude. The amplitude for the production of 
three different particles k, m, n, which is denoted as Mkmn, is the following 

«fem;n(Sl2, z) + a„m;fc(Sl3, z) + akn;m{s23, z) , (67) 

where the first term in the rhs gives the amplitude with the last interaction between 
particles k and m, and so on. The interaction of two particles, for example, n and m, can 
happen in different channels. Here I neglect the isospin ( let / = 0), so the two-particle 
momentum describes the type of interaction. Hence, the amplitude a„m.fc can be written: 

anm;k = Pj^nin;k ) (68) 
J 

where ^^m-fc amplitude with the last interaction of particles n and m with the 

momentum J and isospin / = 0. The integral equation for A^^.j, can be written as 
a generalization of eqs.(^), ( pQ]) and (|52|). I start with the term, where only particles 
m and n interact with each other and the particle k is spectator. The resonance and 
non-resonance production of these particles from the initial state is 

A^".;.(.) = E + (69) 

where C^"^ is the production constant of resonance a and particle k, ^mn-k corresponds to 
the direct production of particles m and n with relative momentum J. The rescattering 
of particles m and n gives the following amplitude: 

E A°/^,J(/-S)-^ , (70) 

(n'm') ^ n'm';nm 
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where the summation over intermediate states with the production of m' and n' particles 
is performed. The expression for triangle diagram can be written in the same way. Finally, 
I get the following three-body dispersion relation: 

{n'm') ^ >n'm';nm („'m') (n"m") ^ J n"m";nm 

/■ ^m'"''('^'l2)^m'^n';m"n"('^'l2) ■S12) . , , , n\ , / / / n\ n n 

J n S12- S12 

where ^ ^ 

N'±';mnis[„ Su) = Af,m'n';mn + E %|^A(.') (72) 

and A(s') is the cutoff function obtained in eq. ( p^ ) with the help of (t)rn'n'{.s). 



4 Dispersion equations for the coupled decay chan- 
nels pp (0~+) Ti^Ti^Ti^, r/Tv^TT^, r]r]7T^, KKtt^ 

The proton-antiproton annihilation in the state J^'" = ^ can originate from the two 
isospin states 1 = and 1=1. Due to isospin conservation in strong interactions there 
exist two integral equations for the following decay channels: 
pp {IJP^ = 10-+) ^ ttOttOttO, ririTT^, KKn^ 
and 

pp {IJ^^ = 00-+) ^ t/ttOttO, KKtt^ 

Let us write integral equations for the decay from 1=1 state. 



a) Reaction pp (10 +j tt'^tt'^tt' 



For this reaction the only difference from the analysis given above is that one should 
take into account the isospin structure of the amplitude. The same analysis has been done 
in refs. [|l^, ||lT|, where the unitarity condition in the irn channel was used to derive the 
amplitude for the rj 37r° decay. This decay goes with violation of isospin symmetry, so 
in both cases initial states have the same quantum numbers. Because of that, here the 
results of refs. [10| and [11] are reproduced, and the contribution from the rjrj and KK 



states is also taken into account. The annihilation amplitude into Stt is 



where 



^2 



(73) 



(74) 



and is the amplitude with the last pions interacting in the isospin state /. It should 

be noted that, due to the C-invariance, there is no term in eq. (|74D with pion interactions 
in the state 1=1. For the channel with isospin / = S- and D-wave interactions should 
be taken into account; this allows one to calculate properly the production of /o and /2 



„0 —400 I p 402 



(75) 
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In the channel with isospin 1 = 2 only S-wave interaction should be accounted for, there- 
fore 

Integral equations for A^'^.^o have the following form: 



TTTT'.TTTT 



TT Si2 - Si2 



2 20 4 

X i{^0.L;A'h3^ + (yaL;^o(s23, z)) j + (-a^^.^o(43, + \v + ^KK, (77) 



where the contributions from the rjr) and KK intermediate states are: 

OJ 

rjrj;7V7r 



(^-^) ^ / } ■ (2 a^ox.^(si3, 



^)). • (79) 



The amphtude a^^.^o in eq. ( |77D can be found only from the annihilation into charged 
pions, so in our approach it could be simply replaced by the direct production amplitude 
of p^TT~ . In the integral equation for A^^.^o only tttt intermediate states are taken into 
account, so one has: 

AZ,.o = \ZM - +{(/-5)-^r X 



X / -, ■ ((a^^;^o(si3, z))q + (-a o(s23, z))q + (a o(s23, 2;))o)(80) 

b) Reactions pp (0~~^) — * ij'qTv^ 



In the ?77r'^ channel the S- and D-wave interaction with the production of ao and 02 
resonances are taken into consideration. The annihilation amplitude is the following: 

M^riwO = ar?r);7r(>('5l2, 2^) + «r?7rO;r?(Sl3, 2;) + a^^o.^(s23, 2;), (81) 
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where 

„ „ _ 400 _|_ p 402 
The integral equation for the A5|;^.^o amphtude has the following form: 



r ^ ^ ^ ^ ^ oj 



X / — -, ^ ((o«7r7r;7ro(Sl3,^))j+ (-5-«7r7r;7ro(«23>^))j+ 

J TT 6]^2 — ''12 "5 y 



4 



where contributions from the 7777 and KK intermediate states are: 

OJ 



A™ = A;;;y|(/-B)-'j + 



I J »7??;»7?? 512 ~ '^12 



OJ 



KK-rm 



i(f R)-A'' f '^''^^ PKK{s'^2)Nl'K;rn^{s'^2, ^12) 
I ) KK:rm J ^lo — S12 



;r;?j TT S12 ~ '^12 

The integral equation for the A^;[o.^ amplitude is: 

ij 



U^-^) f — T-^^i ((S^;7ro(Si3,^))j+ K.rO;^(s23, 

I- J J?7r0;?77r0 J " — S12 

b) Reactions pp (t)~+j — > KKn'^ 
These annihilation amplitudes are: 
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As before, in the KK channel the S- and D-wave interactions and K* resonances in 
Kti channel are accounted for. So, one has: 

For these amplitudes one can get the following integral equations: 



X 



X / -, ■ ((a^^0;K(Sl3>^))^+ («XK;7ro(S23>^))^) ■ 

J TT S10 — Sio 



1 /91 

(is'l2 P7r0if(s'l2)^^0 

■^12 ~ ■^12 

The integral equation for the A^^^^.^o amplitude has the following form: 

OJ r ^ OJ 



X 
KK 



^ f ds[^ Pnn{s'u)Ki,KKis'l2, ^12) ,,2 q U ^ / 2 / / x\ ^ 

J TT S19 — Sio O ' y ' 



■^12 ~ '^12 



4 

where contributions from the rjij and KK intermediate states are as follows: 



I ) 'n'n-KK 



I J ??r);J!'J!: ^ TT •S12 ~ ^12 



L J KK;KK 

f.f m-ll"-" /■ Cig;2 Pgx(g;2)^i^i^;^i^(g;2. ^12) ,^ , 

I ) KK;KK J TT S12 — S12 

A set of integral equations for the reactions pp {IJ^'-^ — 00""^) — > 777r°7r°, KKtt'^ 
may be written in the same way. 
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5 Conclusion 



To summarize, the two-particle discontinuity of the decay amphtude is written, and the 
integral equation is obtained which takes into account not only the S-wave binary in- 
teractions but also interactions with higher angular momenta (P- and D-wave). These 
equations are generalized for the resonance and non-resonance type of two-particle inter- 
actions. A complete three-body dispersion equation is derived, which may be used for the 
analysis of various three-particle reactions. I demonstrate how the coupled three-body 
equations can be written for the 7r°7r°7r°, t/tt^tt^, rjrjir^, KKn^ channels in the pp 
annihilation at rest. 
This work is supported by the grants INTAS-93-0283 and RFFI 96-02-17934. 
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